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ABSTRACT 

We have made a comprehensive transit search for exoplanets down to ~1.5 - 2 Earth radii in the HD 
189733 system, based on 21-days of nearly uninterrupted broadband optical photometry obtained with 
the MOST {Microvariability & Oscillations of STars) satelhte in 2006. We have searched these data 
for reaUstic Umb-darkened transits from exoplanets other than the known hot Jupiter, HD 189733b, 
with periods ranging from about 0.4 days to one week. Monte Carlo statistical tests of the data with 
synthetic transits inserted into the data-set allow us to rule out additional close-in exoplanets with 
sizes ranging from about 0.15 - 0.31 Rj (Jupiter radii), or 1.7 - 3.5 (Earth radii) on orbits whose 
planes are near that of HD 189733b. These null results constrain theories that invoke lower-mass hot 
Super-Earth and hot Neptune planets in orbits similar to HD 189733b due to the inward migration 
of this hot Jupiter. This work also illustrates the feasibility of discovering smaller transiting planets 
around chromospherically active stars. 

Subject headings: planetary systems - stars: individual: HD 189733 - methods: data analysis - 
techniques: photometric 



1. INTRODUCTION 

With the launch o f the MOST (Microvariability & Os 
dilations of STars) (Walker et al.' 2003: Matthews ct al 
[20041 and COROT (Baghn , J003; Barge ct al. 200 " 
satellites, and the upcoming laun ch of the Kepler 
(jBorucki et al.ll2004l: [Basri et al.ll2005l ) space mission, the 
age of transit searches with continuous space-based pho- 
tometry is now upon us. With these new telescopes we 
should now be able to utilize the transit method to probe 
the terrestrial and giant-terrestrial regime of transiting 
exoplanets. We can serious ly anticipate the dis covery 
of a transiting Super-Earth ((Valencia et al.l l2006h sized 
planet orbiting a Sun-like star in the near-term future. 
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Numerical simulations of migrating hot Jupiters 
have indicated that one might expect Earth, Super- 
Earth and Neptune sized planets in similar orbits, 
or in mean-moti o n resonances with the hot Jupiter 
JZhou ct al.' '2005'; 'Ravmond_e£an 120061: iMandell et all 
|2007; Fogg & Nelson 2007). The interior regions of 
these hot Jupiter exoplanetary systems, and specifically 
the mean-motion resonances of these hot Jupiters, have 
been probed with the transit-timing method in a hand- 
ful of transiting exoplanetary systems to date - these are 
TrES-1 (KO yHSteffcn & Agol 2005), HD 209458 (GO 
V) (Agol fc Steffen 1120071: iM iUcr-Ricci ct al. 2007a ), and 
HD 189733 (KO V) Jwinn et al...2007; MiUer-Ricci etaTI 
I2007bl : iPont et aLllsubmittedl ). The transit method, al- 
though currently less sensitive in the mean-motion reso- 
nances of the hot Jupiter, is able to probe a continuous 
range of orbital separations in the interior regions of these 
hot- Jupiter systems. The transit method has only been 
used to sensitively probe the interior regions of one hot 
Jupiter exoplanetary system to date; ICroll et"al] (|2007l ) 
reported a null-result in the HD 209458 system from the 
first space-based transit-search using nearly continuous 
photo metry returned b y the MOST satellite. In that 
work, ICroll et al] (j2007t ) ruled out hot Super-Earth and 
hot Neptune-sized planets larger than 2-4 Earth radii in 
the HD 209458 system with periods ranging from 0.5 
days to 2 weeks. Here we report a search using MOST 
photometry and similar methods ruling out transiting 
exoplanets with periods from 0.4 days to one week in 
the HD 189733 system with orbital inclinations similar 
to that of the known hot Jupiter exoplanct, HD 189733b. 

The existence of a hot Jupiter planet, HD 189733b, in 
th e star system HD 189 733 (V = 7.67) was first reported 
bv lBouchv eFaH (|2005D . The planet was detected from 
radial velocity measurements in the ELODIE metallicity- 
biased search for tra nsiting hot Jup i ters. Photometry of 
HD 189733 taken bv lBouchv et al'l ()2005D indicated the 
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planet HD 189733b transited the star along our line of 
sight. Follow-up work has determined the stellar and 
planetary characteristics of the de tected hot Jupiter to a 
high d egree of accuracy (iBakos et al. 20 06aHWinn et all 
I2006t IWinn et all l2007t iBaines et al.l 1200% . A num- 
ber of researchers have determined characteristics of the 
planet by measuring the drop in flux during the sec- 
ondary eclipse. Th ese results include t he thermal emis- 
sion of the planet (iDeming et al.ll2006f). one of the firs t 
spectra of an extrasolar planet (|Grillmair et al.l l2007f ) , 
as well as the first measure ment of the day-nigh t con- 
trast of an extrasolar planet (jKnutson et al.l l2007y Re- 
cent evidence indicates that HD 189733 is likely a bi- 
nary system with a mid-M dwarf star, HD 189733B (pe- 
riod 3200 yr), as a c ompanion to the main KO star, 
HD 189733 (|Bakos et al. 2006 b). The putative transit- 
ing planets that are searched for in our manuscript are 
assumed to transit the primary star, HD 189733, as does 
the known planet HD 189733b. Radial-velocity mea- 
surements have already ruled out other planets than the 
known hot Jupiter in this sys tem that are mor e massive 
th an approximately ~32 Mm (jWinn et al.ll2006l as quoted 
in lMiller-Ricci et al. I[2007bh . 

Recently, HD 189733 was discovered to display 
quasiperiodic fiux variations at the ~3% level 
(jWinn et al.l l2007i : Matthews et al. in prepara- 
tion). This r e sult i s in concordance with the results of 
IWright et all (|2004D . which indicated that HD 189733 
was a chromospherically active star. Detecting planets 
via the transit method in chromospherically active stars 
is a relatively new topic. In addition to the original 
discovery of HD 189733b aroimd this a ctive star, and the 
pioneering work o f the TEP network (iDeeg et al.lfl998l : 
iDovle et all l2000f ). iHebb et all ([20071) recently found 
that they could rule-out Jupiter and Neptune-sized plan- 
ets with short periods via the transit method around 
the active star AU Mic. We extend these efforts to the 
Neptune and Super-Earth-sized regime by investigating 
the continuous and accurate photometry returned by 
MOST to search for Super-Earth and Neptune sized 
planets with short periods around this active star. 

MOSTs 21-day observations of HD 189733 are pre- 
sented in Matthews et al. (in preparation) . These obser- 
vations have already been used to place limits on other 
hot Super-Earth and hot Neptune exoplanets near the 
mean-motion resonances of the hot J upiter in the system 
through the transit-timing method (jMiller-Ricci et aT71 
[2007b.) . We extend these efforts here by ruling out hot 
Super-Earth and hot Neptune-sized planets using the 
transit method. The rotational modulation displayed 
in MOSTs observations of this system has also been fit 
with a starspot model bv lCroU et al.l ([submittedl ). Their 
model argues in favour of moderate spin-orbit misalign- 
ment between the stellar spin axis and the orbit-normal 
of the hot Jupiter in this system. 

In gl the MOST photometry of HD 189733 is de- 
scribed. The transit search technique is briefly summa- 
rized in Sj3| The Monte Carlo statistics used to estimate 
the sensitivity of the transit search are specified in [JH 
The transit search routine is applied to the MOST HD 
189733 data set in SJU where we are able to rule-out hot 
Super-Earth and hot-Neptune sized planets in this sys- 
tem. Our results are summarized in ^ including a dis- 
cussion of the approximate mass of exoplanets that have 



been ruled out by our results for various compositions of 
extrasolar planets. 

2. MOST PHOTOMETRY OF HD 189733 

The MOST satellite was launc hed on 30 Jun e 2003 and 
its initial mission is des cribed bv lWalker et al.l |2003) and 
iMatthews et"all ([200l . A 15/17.3 -cm Rumak-Maksutov 
telescope feeds two CCDs, one originally dedicated to 
tracking and the other to science, through a single, cus- 
tom, broadband filter (350 - 700 nm). MOST was placed 
in an 820-km altitude circular Sun-synchronous polar or- 
bit with a period of 101.413 minutes. From this vantage 
point, MOST can monitor stars in a continuous viewing 
zone (covering a declination range -|-36° < 5 < —18°) 
within which stars can be monitored without interrup- 
tion for up to 8 weeks. Photometry of very bright stars 
(V < 6) is obtained in Fabry Imaging mode, in which a 
Fabry microlens projects an extended image of the tele- 
scope pupil illuminated by the target starlight to achieve 
the highest precision (see Matthews et al. 2004). Fainter 
stars (down to about V ~ 12) can be observed in Di- 
rect Imaging mode, where defocused images of stars are 
monitored in Science CCD subrasters (see iRowe et al.l 
;2006a). 

MOST observed HD 189733 in Direct Imaging mode 
for 21 days during 31 July - 21 August 2006. For the 
first 14 days of this run, HD 189733 was the exclusive 
target; for the final 7 days, it was shared with another 
Primary Science Target field during each 101.4-min orbit 
of the satellite. Consequently, during the last third of the 
run, HD 189733 was observed only during intervals of 
high stray Earthshine, resulting in somewhat increased 
photometric scatter. The duty cycles of the light curve 
are 94% during the first 14 days and 46% during the 
last 7 days. Early in 2006 the tracking CCD system 
of MOST failed due to a particle hit. Thereafter both 
science and tracking were carried out with the Science 
CCD system. To avoid introducing significant tracking 
errors exposures were limited to 1.5 sec. 14 consecutive 
exposures were then stacked on board to improve signal- 
to-noise, and downloaded from the satellite every 21 sec. 
Approximately 112,000 of these combined observations 
were taken. 

The photometric reduction was performed by one of 
us (JFR) us i ng t echniques similar to those described by 
I Rowe et all ([200 6a) and Rowe et al. (in preparation). 
The reduction procedure of aperture photometry is sim- 
ilar to that applied to groundbased CCD photometry, 
but is non-differential. Our reduction pipeline corrects 
for cosmic ray hits (more frequent during satellite pas- 
sages through the South Atlantic Anomaly [SAA]), the 
varying background due to scattered Earthshine modu- 
lated at the satellite orbital period, and flatfielding ef- 
fects. Further details are given by Matthews et al. (in 
preparation). It should also be noted that HD 189733, 
and its M-dwarf c ompanion are sufficie ntly well separated 
(216 AU or 11.2" [Bakos et al.[ l2006b]). that they do not 
contaminate the photometric signal obtained by MOST. 

2.1. Additional selection and filtering of the data 

For this transit search, additional selection and filter- 
ing of the photometry was done to optimize the data for 
the application of Monte Carlo statistics. The filtering 
is an automated part of the transit search routine, and 
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is thus briefly summarized in ^ This filtering step is 
described in detail here. 

Data obtained during passages through the SAA are 
conservatively excised from the light curve, due to the in- 
creased photometric scatter during those MOST orbital 
phases, without reducing seriously the phase coverage of 
the exoplanetary periods searched. The transits of the 
known giant planet HD 189733b are remove d, at the or- 
bital p eriod P ^ 2.21857d determined by iBakos et al.l 
(|2006af) . The first r^O.Sd of the data-set was excluded 
from the analysis due to an obvious temperature swing- 
in effect common to MOST. The obvious modulation of 
the light curve, du e to starspqt acti vity, was removed 
with a cubic sphne (jPress et al.lll992D . The cubic spline 
was generated with data binned every ^-^500 minutes; this 
trend was then subtracted from the unbinned data. This 
value of '^500 minutes is believed to be optimum, as it 
effectively balances the competing desires of having a suf- 
ficiently short binning time to effectively remove long pe- 
riod modulation, with the desire of having a long binning 
time to ensure that one is able to recover transits with 
periods in the upper range of our interesting period range 
{P > 5d). Limited statistical tests have confirmed that 
this binning value optimizes our sensitivity to transits of 
HD 189733 with periods in the upper range of our period 
range. The signal to noise that an individual transit is 
recovered with can be highly dependent on this binning 
time for transits of relatively long periods (P > 5d). At 
these longer periods the cubic-spline can filter out much 
of the amplitude of the transit signal a low percentage of 
the time. This effect can also be seen with the reduced 
efficiency at longer periods for the ninety-degree incli- 
nation angle case of the Monte Carlo statistics (Figure 

El). 

MOSTsuffers from stray Earthshine with its 101.4-min 
orbital period, and this stray light background can also 
be modulated at a period of 1 day and its first harmonic 
(due to the Sun- synchronous nature of the MOST satel- 
lite orbit). For this reason, sinusoidal fits with periods 
within 1% of 1 and 0.5 d were subtracted from the data. 
After these cuts, any remaining outliers greater than 6a 
were excised. In general, there were few points removed 
by this sigma-cut as few points were such extreme out- 
liers. It is also important to note that the magnitudes 
of the injected transits (Q were always at a level much 
less than this sigma-cut. The data were also median - 
subtracted, for reasons outlined in ICroll et"al] (|2007l ). 
The automatic filtering step removed ~8% of the original 
data (mostly due to the SAA and stray light corrections) . 

The HD 189733 light curves before and after these 
reduction steps are plotted in Figure [1] These fil- 
tered data (and the original reduced data) can be 
downloaded from the MOST Public Data Archive at 
www.astro.ubc.ca/MOST. 

3. TRANSIT SEARCH ALGORITHM 

The transit s earch routine used here is fully described 
in ICroll et"aI1 (|2007t l. and is briefly summarized here. 
The minor modifications to the routine are highlighted 
below. The transit routine consists of five steps: (i) an 
automated filtering step (described in detail in H2.1^ . (ii) 
a modified version of the EEBLS (Ed ge Effec t Box-fitting 
Least Squares) algorithm of Kovacs et al.l ([2002 ) that 
searches for box-shaped transits, (iii) selection criteria 
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Fig. 1.— The 2006 (A) HD 189733 MOST data-set before the 
automated filtering step. The data following (B) the automatic 
filtering step, in which the transits of the known planet were re- 
moved, and a cubic spline removed the obvious long-term variabil- 
ity. (C) The same data are shown binned in 30 minute intervals at 
a different vertical scale. 



to select the the best candidates that can be tested for 
astrophysical plausibility, (iv) realistic transit-fitting and 
refinement of the transit parameters, and finally (v) de- 
tection criteria to differentiate false positives from likely 
transit candidates. 

The detection criteria used in this study are summa- 
rized here. We use the transit to anti-transit (the best-fit 
brightening event with the shap e of a transit) ratio statis- 
tic (RS) of iBurke et al.l ()2006[ ) to quantify a behevable 
transit. The RS value that a transit has to exceed to be 
deemed a believable transit is: Ax^%/Ax?.% > 0.735, 
where Ax^ and Ax^ refer to the normal (the sum 
of sqaured residuals) for the best-fit transit, and anti- 
transit, respectively. This value was motivated by the 
results of the Monte Carlo tests (12]). The RS statistic 
serves to objectively quantify the significance of the mea- 
sured signal against the correlated scatter in the data, 
and to thus stringently rule-out false positives. The cor - 
related signal to noise calculation of iPont et al.l ()2006f ) 
serves a similar purpose. The RS will be used in this 
present application. 

The restriction on the quantity f, a measure of the 
ratio of the strength of any one s pecific trans it com- 
par ed to the others, a s described in ICroll et al.l ([2007), 
an d IBurke et alllpOOl . was also relaxed such that tran- 
sits were accepted as long as f < 0.85. This motivation 
for this less stringent criterion was based on the Monte 
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Carlo statistics that indicated that a more conservative 
criterion on f would rule out realistic transits, without 
an appreciable decrease in the number of false positives 
in our data-set. 

The precise values used for various parameters in the 
EEBLS step (ii) in the transit search routine are sum- 
marized in Table [TJ We have decreased the minimum 
value probed by the EEBLS routine in the application 
of this transit search routine to HD 189733, compare d 
to its application to HD 20945 8 in ICroU ct al. (' 20071 ). 
This is due to the conclusions of lSahu et al . (2006)), who 
observed a number of ultra-short period (P < 1.0 d) 
planet candidates in the SWEEPs survey, and noted that 
these planets may preferrentially be found around low- 
mass stars. HP 189733 has a stellar mass of 0.82 Mq 
(iBakos et all [2006a). below the hmit of 0.88 Mq sug- 
gested bv ISahu et al.l ((2006) to be indicative of a pref- 
erence for ulta-short period planets. The minimum pe- 
riod investigated by the EEBLS routine is thus now OAd, 
a value corresponding to a semi-major axis, Op, of just 
under three times the stellar radius (op « 2.5i?*). For 
comparison, OGLE-TR-56b has ap«4.4i?*, the smallest 
orbital radius in terms of its stellar radius of a confirmed 
planet to date. The maximum period of 6.9d is cho- 
sen similar to th e maximum value for HD 209458 in 
iCroll et all (l2007h . as a value that would allow one to 
observe three distinct transits in the 21 d data-set. This 
period range of 0.4 to 6.9 d corresponds to 0.01 to 0.07 
AU. The minimum fractional transit length is 0.75 times 
that of the calculated average fractional transit length 
at that period - a value chosen to ensure the routine re- 
mains adequately sensitive to transiting planets with non 
edge-on inclination angles. We also increase the number 
of bins, Nb, used by the EEBLS algo rithm for short pe - 
riods compared to the application of ICroll et al.l ()2007f ) 
to HD 209458. The number of bins, Nb, used by the 
EEBLS algorithm for each period, P, is determined by 
the following formula: [0.7 x exp(— P) -I- 1.0]x20.0/Qmi, 
where Qmi is the minimum fractional transit length to be 
tested. Qmi, and the maximum fractional transit length, 
Qma, are determined as given in Table (TJ from our av- 
erage estimate of the fractional transit leng th for a hy- 
pothe tical planet, Qmp, at various periods (jCroll et al.l 
[20?)1 . 

In our transit selection step we now select the top 
three EEBLS SR signals, instead of the top two, as 
well as up to two transit signals with SR signals that 
are at least 1.8 times the noise flo or. A full explana - 
tion of these details is provided in ICroll et all (|2007D . 
The stellar characteristics used to determine the realistic 
limb-darkened transit model (iv) are summarized in Ta- 
ble m The precise transit parameters are refined using 
a Marquardt-Lev enberg routine (described in detail in 
ICroll et al.ll2007t l. The method used to refine our tran- 
sit parameters is very similar to the Newton-R aphson 
method apphed bv lColher-Cameron et alj {j20Q6D for the 
same purpose. 

4. MONTE CARLO STATISTICS 

4.1. Monte Carlo setup 

To assess the sensitivity of our routine and the MOST 
data-set to smaller transiting planets in the HD 189733 
system, simulated transits for planets of various radii and 
orbital parameters were inserted into the MOS'Tphotom- 
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Fig. 2.— The EEBLS spectrum of the MOST HD 189733 data- 
set, plotted in frequency (top), and period (bottom). The 3 can- 
didates that passed the transit selection criteria arc shown by the 
solid vertical lines. 

etry and Monte Carlo statistics of the transit recovery 
rate were generated. We thus proceed under the assump- 
tion that all signals in the data-set are noise, and we place 
as sensitive of limits as possible under this paradigm. Re- 
alistic limb-darkened transits due to planets with various 
radii, Rp, orbital phases, (j), periods, P, and inclinations, 
i, were inserted into the 20( 36 MOST HD 189733 data. 
The full non-linear model of iMandel fc Agoli (|2002D was 
used to ensure the inserted transits were as realistic as 
possible. Circular orbits were used. The M-dwarf com- 
panion star, HD 189733B, would have a negligible effect 
on the orbit of a close-in putative planet in orbit around 
HD 189733, and thus the effect of this binary companion 
can be safely ignored. 

These modified data were then subjected to the transit 
search algorithm as discussed in [J3l It should be noted 
that this includes the automated filtering step, and thus 
the synthetic transits were injected before the same filter- 
ing process our original data was subjected to, therefore 
ensuring the validity of the following Monte Carlo limits. 
Transits were in serted with logarithmic period spacing 
(as discussed in ICroll et al.l [2007D . with r]inp = 0.04 in 
the period range 0.45d < Pinp < 6.7d. In total, 36 pe- 
riod steps were used for the 90° inclination angle cases. 
For the 86° and 82° inclination angle cases transits were 
inserted with this same logarithmic period spacing until 
the period exceeded the maximum period, or semi-major 
axis, that would produce a transit, a < (i?* + Rp)/cosi. 
For each trial period, simulated transits corresponding to 
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TABLE 1 

EEBLS INPUT PARAMETERS 



definition 



value 



Np Number of Period Points searched 

77 Logarithmic period step 

P-min Minimum period submitted to EEBLS algorithm 

Pmax Maximum period submitted to EEBLS algorithm 

Qmi Minimum fractional transit length to be tested 

Qma Maximum fractional transit length to be tested 

Nb Number of bins in the folded time series at each test period 



15 000 
0.000219 
OAd 
6.M 

O.lbxQmp 
l.lxQmp 
(0.7 X exp(-P) + 1.0)x20.0/Qmi 



TABLE 2 

System parameters 



var 


definition 


value 




Stellar radius 


0.755 Rq 


M. 


Stellar Mass 


0.82 Mq 


Cl 


Non-linear limb-darkening parameter 1 


0.82523168'^ 


C2 


Non-linear limb-darkening parameter 2 


-0.97755590^^ 


C3 


Non-linear limb-darkening parameter 3 


1.6867374<: 


C4 


Non-linear limb-darkening parameter 4 


-0.65702013'^ 



^ Parameter obtained from IPont et a T ( submitted')'' Parameter 
obtained from lBakos et ah I I l2006al)'= These limb-darkening param- 
eters are fully described in IMandel &: Agoll (|2002f l and obtained 
from lMiller-Ricci et al"~l Il2007bl) . 



20 different exoplanet radii were inserted, sampling the 
period-radius space of interest. For each period and ra- 
dius at least 100 phases were inserted. For each of these 
points, the phase </) was generated randomly to be in the 
range < < 1. Because a 2.4 Ghz Pentium processor 
with 1 Gbyte of memory can perform the transit search 
algorithm on an individual MOST HD 189733 data set in 
^1 minute, exploration of the entire grid just mentioned 
for all inclinations involves ~3xl0'^ iterations and 0.6 
CPU years. The calculation was performed on the LeV- 
errier Beowulf cluster in the Department of Physics and 
Astronomy at the University of British Columbia using 
46 dual-CPU compute nodes. 

An inserted transit was judged to be detected if the 
parameters (/> and P returned by the transit search algo- 
rithm were sufficiently close to the input values: 0i„p and 
Pi„p. The returned period had to satisfy the following 
criteria: \P - P„p| < O.OSd and |P/P™p - 1| < 1%. The 
limits of the criterion for were dependent on the plane- 
tary orbital period, because as the period of the putative 
transiting planet decreases, the fractional transit length 
increases accordingly. Thus the accuracy required in the 
determination of was relaxed for shorter periods. The 
criterion on is: |</) — </)mp| < 0.09 — 0.0106 x (Pi„p - 
0.4d)/d. This corresponds to accepting an error in phase 
of up to 9% at 0.4c;, but only up to 2% at Id. Obvi- 
ous multiples of the period of the inserted planet, up to 
4 times the inserted period, Pi„p, as well as half-period 
(P « jPmp), and one-third period (P « \Pin-p) solutions 
were also accepted. For the Monte Carlo tests when a 
transit was deemed recovered the correct period was re- 
covered 96% of the time, while the remainder returned 
a period that was double, triple, quadruple, one-half, or 
one-third of the inserted period, Pmp. 

The level above which a transit recovered from the 



data-set could be considered believable was also deter- 
mined due to these Monte Carlo tests, at a level of 
Ax^%/Ax?.% > 0.735 that ruled out 95% of the spuri- 
ous candidates. Spurious candidates are defined as those 
candidates where the parameters returned by the transit 
search routine are different from the parameters inserted 
for the Monte Carlo tests. The motivation for this cri- 
terion is to ensure that a systematic event is misidenti- 
fied as a believable transiting planet less than 5% of the 
time. Interestingly we can rigorously rule out these spu- 
rious candidates, with our criterion set at a level where 
the most significant transit candidate can be less signifi- 
cant than the most significant anti-transit (brightening) 
event. Fur t her d iscussion on this criterion is given in 
ICroll et "all ()2007f ). So as to not unduly skew this statis- 
tic we only generate this statistic using putative planets 
with radii large enough that the routine has a realis- 
tic chance of detecting these planets. We judge this to 
be planets greater than approximately the 25% contour 
of our Monte Carlo statistics without the RS criterion. 
Thus only planets with radii greater than the following 
limit are used: Pp„ip > (0.10 + ^ x (P„ - 0.4d)/d) 

4.2. Monte Carlo results 

The fractions of times that the artificially inserted 
transits were recovered from the data for various radii, 
periods, and inclination angles are given in Figure[3) Also 
indicated in Figure[3]is the 68% contour limit that would 
be placed without the Ax^%/Ax?.% > 0.735 criterion. 
The close agreement between the 68% contour with or 
without this criterion indicates that this criterion does 
not seriously impact the sensitivity of the Monte Carlo 
statistics generated, while providing a robust limit so as 
to avoid false positives and spurious detections. 

We also explore the distribution of spurious signals 
across period-radii space. These values are displayed in 
Figure HI and indicate that spurious signals are negligi- 
ble (below 1%) other than the intermediate areas where 
the routine has a small to moderate chance to correctly 
determine the inserted transi t. This resu l t is e xpected, 
as it is similar to the result of ICroll et all (|2007f ) for HD 
209458. 

The performance of our routine is also expected to be 
slightly degraded at harmonics and subharmonics of the 
periods that the sinusoids were removed from the data 
at (P=1.0d, and P=0.5d). Due to the fact that the frac- 
tional length of a transit is comparitively large for short 
periods of 0.5 - 1.0 d, it can be difficult to completely 
differentiate and properly recover a transit signal near 
the harmonics and subharmonics of these periods in all 
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cases. The period ranges deraonstrating these drops in 
survey sensitivity are very narrow ~ 2%). Limited 
statistical tests confirm that the loss in sensitivity near 
these periods is qu ite small, similar to that of HD 209458 
(jCroll et al.ll2007D . 

The harmonics and subharmonics of the orbital pe- 
riod of the known planet HD 189733b, P « 2.218573d 
(|Bakos et al.ll2006al ). are also expected to show decreased 
sensitivity to transits due to the fact that the transits 
of the known planet were excised from the data. This 
reduction of sensitivity at these periods is unfortunate 
considering that it is expected that low-mass planets 
show a preference f or these mean-motion resonant or bits 
Jriiommcs 2005; Papaloizou fc Szuszkiewicz1l2005l and 
IZhou et al.ll2005( ). The fractional length in phase that 
was removed at P ~ 2.218573 to exclude the transit of 
the known planet was 4.0% of the total. Thus for the 
subharmonics of the known planet (P ~ 4.44d, 6.66d) in 
approximately 4.0% of the cases the routine will be un- 
able to recover a putative transiting planet, as the tran- 
sits will coincide with the transits of the known planet, 
and this information will have been completely removed. 
For the harmonics of the orbital period of HD 189733b, 
the situation is better, since not all of the transit events 
would have been removed from the data. For periods 
near the first harmonic (P « 1.1 Id), every second occu- 
rance of a transit will be lost in only 8.0% of the cases. 
For the second harmonic (P « 0.74(i), every third tran- 
sit would be missed in 12.0% of the cases. For periods 
near these values the sensitivity limits can be expected 
to be only slightly degraded from those shown in Fig- 
ure [3l Even for exact harmonics the routine's sensitivity 
should be only marginally worse than the limits quoted in 
Figure [3l as the remaining transits should allow the cor- 
rect period (or a multiple of the period) to be recovered. 
The mean-motion resonances from the known exoplanet 
HD 189733b have been examined using the tr ansit-timing 
techni que with ground-base d photomet r y by I Winn et al.l 
(|2007l ). HST photometry bv lPont et all (Isubmittedl). an d 
with MO^T photometry bv iMiller-Ricci et al. I (|2007bf ). 
It should also be noted that the detection limits pre- 
sented here are valid for circular orbits, but are largely 
applicable to orbits of other eccentricities, as summa- 
rized for the analogous case of HD 209458 in lCroll et all 
(|200l . 

We have made the simplifying assumption that the hy- 
pothetical planet that transits the star to produce the 
synthetic transits of our Monte Carlo statistics does not 
pass over prominent starspots along our line of sight. The 
known planet has not been found to occult a prominent 
starsp ot in the MOST photometry (|Miller-Ricci et al71 
l2007bf ). although a preliminary investigation of HST pho- 
tometry of HD 189733 indicates that the planet does oc- 
cult small to moderate sized star spots during their ob- 
servations (jPont et aLllsubmittedl ). Given that the rota- 
tional period of the star is greater than the period range 
that is investigated here, the limits presented here would 
likely only be slightly adversely affected if this hypothet- 
ical planet does pass in front of a prominent starspot. 
This is because the impact on the transit dip of the planet 
occulting a large starspot would be relatively short, and 
should thus not severely affect our limits. 

Thus the search routine described above in Sj3] should 
be able to detect planets with radii greater than the lim- 



its given in Figure [3] with the MOST data. In the most 
optimistic case of an edge-on 90° inclination angle tran- 
sit, for planets with periods from approximately half- 
a-day to one week, this limit is approximately 0.15 to 
0.31 Rj, or 1.7 to 3.5 P®, respectively, with 95% confi- 
dence. These limits would scale up, or down linearly if 
the expected stellar radius, P* = 0.755 Rq, is an over- or 
under-estimate, and are thus perhaps better expressed as 
Pp/P* = 0.020 and Pp/P* = 0.042 for periods from half- 
a-day to one week, respectively. For periods at the upper 
range of those investigated here (P > 6.5d) the routine is 
unable to pick out the transit in a low percentage of cases, 
because the cubic spline has removed significant power at 
that period; this effect can be observed as a reduction of 
sensitivity for longer periods in the ninety-degree inclina- 
tion angle {i—9Q°) panel of Fig. \3\ The performance of 
our automatic routine in this period-radius space prob- 
ably underestimates the true sensitivity of our survey, 
as it would be obvious that the cubic spline has over- 
corrected in these cases through an individual inspection 
of the phased curve. Also, limited statistical tests have 
shown that marginally improved limits could probably 
be set for the shortest periods (P < 3d) by excluding the 
last seven days of observations with increased noise, and 
only producing these limits with the first 14 days of the 
highest quality data. 

5. HD 189733 TRANSIT SEARCH 

The MOST HD 189733 2006 data-set was submitted 
to the analysis outlined in fj3l There were no transiting 
planet candidates that met the detection criteria. As the 
Monte Carlo tests indicate we would be able to detect 
exoplanets in this system with periods from half-a-day 
to one week with radii larger than the radius limits given 
in Figure [21 we can safely rule out these sized-exoplanets 
in this system. 

We also present the details of the candidate with the 
greatest improvement in even though evidence for 
its existence is marginal. The phased diagram of this 
candidate (Figure [5l) is not particularly convincing, and 
this candidate did not meet the RS detection criterion, 
{Ax^%/Axl% « 0.438 ^ 0.735). For these reasons we 
do not report it as a probable transiting candidate. The 
specifics of this putative candidate are given in Table [H 
This event, with a period of P « 0.752d, has modest 
statistical significance, and thus we report it is a possi- 
ble, but unlikely, candidate. This candidate would cor- 
respond to a putative hot Super-Earth, with period P = 
0.752(i, radius Rp = 0.134Pj (1.50P©), and incHnation i 
— 86.0°. Interestingly the period of this planet places it 
quite close to the 1:3 interior mean- motion resonance to 
the known planet (P--0.74d). For p « 3000 kg m"^ the 
mass of this putative planet would be approximately 1.8 
M® (0.006 Mj). A transit time combined with the sup- 
posed period is given in Table[3l Evidence for a transit at 
this period is marginal, but additional MOST photome- 
try of the HD 189733 system should confirm or deny its 
existence. 

The most significant brightening event was one ob- 
served with a period of approximately P « 6.28d. It 
is likely statistical in nature, and thus is not expected to 
be related to any specific astrophysical process. 

Although our transit search covers a continuous set 
of periods, we do not expect the existence of another 
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Fig. 3. — Confidence Umits of transit detection as a function 
of planet radius and orbital period, for different orbital inclina- 
tions, based on Monte Carlo statistics. The crosses represent the 
radii and periods at which synthetic transits were inserted in the 
data. The dotted, thick-solid, and thin-solid lines represent the 
99, 95 and 68% confidence contours, respectively. The thick dot- 
dashed curve represents the 68% confidence contour if the criterion 
Ax^%/Ax?.% > 0.735 was not used. The near vertical dot-dash 
line in the later panels indicates the maximum period that produces 
a transit for that given inclination angle. Note the logarithmic pe- 
riod scaling on the x-axis. At least 100 phases were inserted for 
each of the radii-phase points. 
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planet with a semimajor axis 'too close' to that of the 
known planet HD189733b, since the gravitational per- 
turbations of the known Jupiter-scale planet will desta- 
bilize the orbit of another object. Taking the known 
planet's data - ab=0 .0313 AU, and mass Mb— 1.15 Mj^p 
(jBouchv et alJ 12005) in orbit around a s tar of mass 
M*=0.82 solar masses (|Bakos et alJ l200^a[) - orbits in 
an annulus ai„ = (1 — A)ab < Of, < (1 -|- A)af, — Uout will 
be unsta ble. The fractio nal width A of the annulus is 
given by (|Gladmanlll99l : A = 2.40(Mf,/M,)i/3 where 
the assumption is being made that the mass of the sec- 
ond planet is much less than the known planet (if not, 
the width of the zone scales as the cube root of the sum 
of the masses of the planets). Here A=0.26, and so we 
calculate that planets with periods in the range 1.41 - 
3.14 days will be intrinsically unstable due to the pres- 
ence of the known planet. Thus, some of the parameter 
space searched in Figure[3]was very unlikely to have been 
occupied by a second planet - Trojans lagging or leading 
(iFord fc Gaudill2006f ) HD 189733b being a notable ex- 
ception. 

We can also exclude transits from planets with longer 
orbital periods. There were no transits from edge-on {i = 
90°) planets with radii greater than ~0.35 Rj (3.9 Rq) 
during MOSTs observations of HD 189733, as a planet 
of this size would be readily visible with even a cursory 
inspection of the data. A planet of this size would cor- 
respond t o ~ 2 Mq, as calculated from the models of 
iFortnev I (|2007t l assuming a composition similar to Nep- 
tune. We can also place a limit on transiting planets with 
orbital periods from 7 d < P < 10.5 d - in this region we 
would expect two transits in our data-set. Limited statis- 
tical tests indicate that we can rule out planets with radii 
greater than approximately ^^0.32 Rj (3.6 Rp) iri this re- 
gion - this corresponds to ^^14 Mq ()Fortnev Il2007f ). again 
assuming a composition similar to Neptune. We do not 
provide a formal limit from Monte Carlo estimates for 
this range of periods as we prefer to restrict our formal 
transit search to periods that will result in three transits 
in our data-set. 

As an additional sanity check this transit search 
method was applied to the current data-set without the 
removal of the transits of the known planet, HD 189733b. 
As the transits are apparent in the unbinned data (Fig- 
ure [T|) , this causes the cubic-spline to over-correct and 
remove a moderate portion of the signal at this period. 
Even so, the routine correctly uncovers the transit of 
HD 189733b with reaso nable accuracy. Th e actual pa- 
rameters of HD 189733b (jBakos et al■ll2006a^ were recov- 
ered within 0.04% and 0.3% in period and phase, respec- 
tively, while the inclination angle and radius were recov- 
ered within 1.0° and 0.20i?,/ of the actual parameters. 
The moderate discrepancy between the inserted and re- 
covered radius is due to the fact that the cubic spline 
overcorrects and removes a portion of the transit-signal. 

6. DISCUSSION 

MOSTs 2006 observations of HD 189733 have been 
searched for evidence of other exoplanets in the sys- 
tem. T he transit-search routine described in ICroll et aTl 
()2007| ) has been used to search for realistic limb-darkened 
transits to take advantage of the precise, near-continuous 
photometry returned by MOST and the fact that the 
stellar characteristics of the star, HD 189733, are well 



TABLE 3 

Transit candidate of marginal significance 



parameter 



value 



Rp 

i 

Ephemeris Minimum (JD-2451545) 
Mass (assuming p ^ 3000 kg m~'^) 



0.7517d 
1.50R(B (0.134_Rj) 
86.0° 
2403.711 
0.438 
1.8 Me (0.006 Mj) 



-0.0100 - 



_ o.oooo 




0.7 0.75 O.J 



-0.001 - 



0.0005 - 



Fig. 5. — The best transit candidate as identified by our analysis 
of MOSTs HD 189733 2006 data-set, but with marginal signifi- 
cance. The thick red-line represents the transit model, while the 
thin green-dashed line represents the constant brightness model. 
At top, the unbinned data are shown in blue, and the binned data 
are shown in black; at bottom the binned data only. As this candi- 
date failed the improvement in transit over anti-transit ratio statis- 
tic (RS) criterion, we report it as a possible, but unlikely candidate. 
The period and radius of this putative planet would be approxi- 
mately 0.752d and 0.134i?,/ (1.50i?e). 

established. Monte Carlo statistics were generated using 
this routine and indicate that this routine in combina- 
tion with the aforementioned MOST data have placed 
robust limits in the size of transiting bodies that have 
been ruled out in this system for a range of periods 
and inclination angles. In the most optimistic case of 
edge-on transits, planets with radii greater than 0.15 Rj 
(~1.7 i?®) to 0.31 Rj (~3.5 Rq) with periods from half- 
a-day to seven days, respectively, have been ruled out 
with 95% confidence through this analysis. For orbits 
co-planar with that of the known planet (i « 86°) the 
limits are similar, except for longer periods {P>4d). For 
these longer periods the planets cease to transit the star 
along our line of sight and thus we cannot place firm lim- 
its for periods greater than 4-days at z ~ 86°, or for i < 
86°. Specifically, we have been able to rule out transit- 
ing planets in this system with radii greater than those 
given in Figure [31 This work has constrained theories 
that predict hot Super-Earths, and hot Ncptunes due 
to the inwar d migration of HD 189733b, such as thos e 
proposed by IZhou et all (I2005D. iRavmond et al.l (120061). 
iMandell et al.l (|2007D . and lFogg k Nelson I (|2007D . 



Looking for Super-Earths in the HD 189733 system 



9 



It should also be noted that we are able to place very 
sensitive limits on the size of planets that can be ruled 
out in this system despite the prominent ~3% rotational 
modulation observed in this system. By removing the 
observed rotational modulation using a cubic spine with 
a binning time longer than the expected duration of the 
transit of the planet across the star, we are able to re- 
move the effects of the prominent modulation while re- 
maining sensitive to the transits of planets as small as 
those expected from hot Super-Earths and hot Neptunes. 
These methods should be applicable to current or future 
space-based transit search missions, such as COROT and 
Kepler. 

We are also able to place limits on the size of Trojan 
planets that have been ruled out in this analysis, assum- 
ing the Trojan consistently transits the star with an in- 
clination angle close to that of the known planet (i~86°). 
We us e the definition of Trojans given bv lFord fc Gaiidil 
(|2006f ) as objects occupying the L4 or L5 stable Lagrange 
points. Trojans lagging or leading HD 189733b with a 
radius above R = 0.24 Rj (2.7i?®) should have been de- 
tected with 95% confidence. Using a mean density of 
jO« 3000 kg m~^ this corresponds to 11 M^. A detailed 
study of this photometry should be able to radically im- 
prove this limit on the size of Trojans that can be ruled 
out with this photometry, assuming the Trojan consis- 
tently transits the star along our line of sight. 

6.1. Mass constraints on other Exoplanets in the system 

The above limits on the smallest planetary radii ex- 
cluded by MOST for transiting planets allow some esti- 
mates on the type of planets that are being excluded in 
this system. For Gas Giant planets, including hot Nep- 
tunes, radii change with time as the planet cools, so we 
make use of the fact that t he age of the HP 189733 sys- 
tem is greater th an 1 Gvr (iMelo et al.ll2006l ). We follow 
the suggestion of lWinn et al.l ( 20071 ) and adopt an age of 
~4.5 Gyr. We thus can refer to cooling models similar 
to the planets in the Solar System. One exception would 
be the consideration of possible planets in orbits smaller 
than that of HD189733b, where tidal heating will lead to 
larger radii per given mass and composition. 

Our radius hmit of 1.7 - 3.5 i?® already places us in the 
realm of Super-Earth and hot Neptune planets. iFortnevI 
(|2P07) has recently provided theoretical predictions for 
exoplanetary radii for a diverse range of possible exoplan- 
etary compositions, including hot Super-Earth , and hot 
Neptu ne planets. The Super-Earth models of iFortnevI 
(|2007f) ag ree well with t h e mo r e det ailed Super-Earth 
models of Valenci a et al.l ()2006l . l2007f) . For our Super- 
Earth exoplan etary limits we use equation 7 and 8 of 
iFortnev I ()2007ll . Our inner period limit oi Rp < I.IR® 
allows us to rule out planets more massive than 1.3M0 
for a pure ice planet, 4.8M0 for a pure rock planet, and 
27M0 for a pure iron planet. Our outer period limit 
of Rp < 3.5i?0 only allows to place a useful limit on 
pure ice Super-Earth planets. At this limit we are able 



to rule out planets more massive than 23M0 for pure 
ice planets. Pure rock and pure iron planets of this size 
are moderately less massive, and comparable in mass, re- 
spectively, to HD 189733b itself - planets of this mass, 
of course, have already been ruled out in this system 
with radi al-velocity measurements ^Winn "erani2006. as 
quoted in iMiller-Ricci et al. II2007H) . We note that our 
limit on ice planets may be more sensitive than the limits 
quoted here as one would expect that ice planets at such 
small orbital periods would have steam atmospheres, and 
would thus be significantly inflated at these small orbital 
seperations (jKuchner Il200l IFortnev II2007D . 

We are also able to place limits on the mass of hot 
Neptune-type planets that have been ruled out in our 
analysis. Hot Neptune planets should have cores simi- 
lar in density to giant ice-plane ts, and an outer atmo- 
sphere of H/He (jFortnev I [20071) - they should thus be 
moderately less dense than giant-terrestrial ice planets, 
and thus be o f a similar den sity to Neptune in our own 
solar system (' Fortnev1l2007D . Our inner period radius 
hmit {R < 1.7i?0) excludes hot Neptunes entirely, as hot 
Neptune exoplanets are expected to be of similar size to 
Neptune (3.9 Ra,) (IFortnev 1 120071: iValencia et alll2006l : 



I Valencia et al.ll2007f ). At our outer period radius limit 
{R < 3.5i?0) we can rule out hot Neptune exoplanets 
more massive than 12.2 M®. 

The radius limits presented here, combin ed with the 
mass - limits from transit-timin g analyses (I Winn et al.l 



I2007t IMiller-Ricci et al. Il2007bl : iPont et aLlfsubmittedl ). 

place very firm constraints on the sizes and masses 
of bodies that could still reside in close orbits to HD 
189733. Re-observations of this system planned with 
MOS'T should be able to further reduce the size of planets 
that could remain undetected in this system, and further 
constrain theories that predict hot Super-Earth, and hot 
Neptune planets in similar orbits to that of the known 
hot Jupiter exoplanct HD 189733b. HD 189733 thus 
joins the gro wing list of h ot Jupiter exoplanetary systems 
(HD 209458 ICroll et all 2007; Miller- Ricci ct al. 2007;^; 
lAgol fc Steffen II2007| . TrES-1 tSteffen fc Agolii2005 ]) for 
which firm constraints have been placed on the size of 
other exoplanets in similar orbits to that of the known 
hot Jupiter by a variety of methods. 
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